Suppressed compressibility at large scale in jammed packings of size disperse spheres 
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We analyze the large scale structure and fluctuations of jammed packings of size disperse spheres, 
produced in a granular experiment as well as numerically. While the structure factor of the packings 
reveals no unusual behavior for small wavevectors, the compressibility displays an anomalous linear 
dependence at low wavectors and vanishes when q — > 0. We show that such behavior occurs because 
jammed packings of size disperse spheres have no bulk fluctuations of the volume fraction and are 
thus hyperuniform, a property not observed experimentally before. Our results apply to arbitrary 
particle size distributions. For continuous distributions, we derive a perturbative expression for the 
compressibility that is accurate for polydispersity up to about 30%. 

PACS numbers: 05.10.-a, 05.20.Jj, 64.70.Pf 



When an assembly of hard particles is compressed, 
there comes a point where further compression is diffi- 
cult because the required pressure is too large. A similar 
"jammed" state can be obtained with soft repulsive par- 
ticles (as in emulsions or foams), at a particular volume 
fraction <j) c above which the stability of the packing is 
controlled by the elasticity of the particles. A large body 
of recent work in both theory and experiment has char- 
acterized the properties of jammed packings (HQ. 

While much attention has focussed on contacts at the 
interparticle scale, and force networks and connectivity 
at larger scales [l|-|3(, there has been comparatively lit- 
tle research into fluctuations and response at very large 
scales. In Ref. Q the low-wavevector behaviour of the 
structure factor, S(q), in a monodispcrsc system of hard 
frictionless spheres was studied numerically, revealing 
"unexpectedly" weak density fluctuations at low q, 



That this is a highly topical question is demonstrated by 



S(q) « aq, 



(1) 



for some constant a, as was also found later for soft par- 
ticles 0. This should be compared to the behaviour in 
liquids (including hard sphere fluids at (f> < 4> c ), where 
S(q) w S(q -> 0) + a'q 2 [ff . These results imply that bulk 
fluctuations in the number density are suppressed at </> c : 
in d dimensions, fluctuations of particle number in sub- 
systems of linear size L scale as (AiV 2 ) L ~ L d ~ x . Such 
suppressed density fluctuations are the defining feature 
of "hyperuniform" materials 

Two recent papers reported surprising results, failing 
to detect the behaviour in Eq. ([1]). A numerical simula- 
tion of binary mixtures @ found that © only holds for 
the particular case of a monodisperse system. Similarly, 
a confocal microscopy study of a jammed system of 
moderately polydisperse colloidal hard spheres also failed 
to observe Eq. ([T]). Both studies report that S(q) in size 
disperse systems is different from Eq. (JlJ, and suggest 
that size disperse packings might not be hyperuniform. 



the recent work of Ref. [Ill ] , where independently of our 
approach hyperuniformity was detected, using two-point 
probability functions. 

Here, we report the first experimental observation of 
vanishing fluctuations of the volume fraction and hyper- 
uniformity in a granular experiment. The same observa- 
tion is made for numerically produced polydisperse pack- 
ings with arbitrary size distributions. In contrast with 
0, [Tol we consider not only the structure factor S(q), 
but also the isothermal compressibility, xt(<z), the latter 
being central to our analysis of size disperse packings. 
We propose a novel perturbative approach for extracting 
Xr(q) for continuous size distributions, which in general 
is a non-trivial task, and explain the connection between 
vanishing compressibility and hyperuniformity. 

First, we briefly describe our systems. Experimentally 
we produce dense random granular packings by slowly 
compressing horizontally vibrated bidisperse brass disks. 
Typically 4500 large disks (diameter 5 ± 0.025 mm) and 
3500 small disks (diameter 4 ±0.025 mm), surrounded by 
rigid walls, arc placed on an oscillating glass plate (am- 
plitude 5mm, frequency 10 Hz). The packing fraction 
is increased logarithmically slowly (d<fi/d\og(t) ~ 10~ 2 ) 
until the force F measured on the compressing wall in- 
creases sharply from F — 0.05 Mg to F > Mg, where 
M = 2 kg is the typical total mass of the "grains" . At 
that point the packing jams, grains stop moving and the 
force measured at the wall remains finite when the vi- 
bration is switched off (see Ref. [l3| for more details). 
We take a high resolution picture of the entire packing 
(2048 x 2048 pixels) at the largest packing fraction, cor- 
recting for optical distortion. We detect the positions of 
the grains with a resolution of 20 /xm and retain ps 4000 
large grains and ~ 3000 small grains at least 1 cm away 
from the walls. We analyze two independent packings, 
produced starting from uncorrelatcd initial conditions. 
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Rattler particles are always included in the analysis. Nu- 
merically we generate 3c? polydispcrsc sphere packings at 
4> c using soft repulsive particles, as first proposed in 
using conjugate gradient methods and small decompres- 
sion steps to prepare packings exactly at 4> c P, [l2| ■ For 
each set of parameters, we prepare a single, very large 
configuration composed of N = 64, 000 particles. We 
study both a 50:50 binary mixture of spheres with di- 
ameter ratio R 6 [1,2], or systems with a continuous size 
distribution, which we take as a flat distribution centered 
around the average value a. We studied polydispcrsitics 
up to p — 0.4, where p is the standard deviation of the 
size distribution divided by a . 

We now provide some definitions. Consider a size 
disperse system composed of n species, containing Ni 
particles of species i, with diameter CTj. We define 
N = Yli=i Ni, the density p = N/V, the concentra- 
tions x.i = Ni/N, and the partial density fields Pi(q) = 
Sj^i elqri ' where Tj is the position of particle j be- 
longing to species i. The partial structure factors read 
Sij(q) = jj {pi(q)pj(—q)), and we collect them in a ma- 
trix, S(q). The total structure factor is defined as usual: 



(2) 



i=l j=l 



When we analyze 5(g) in our experimental granular 
packings, see Fig. [TJt, we find that the low-g behaviour 
is not compatible with Eq. ([T]). The same observation 
holds, see Fig. [^t, for the evolution of 5(g) for numerical 
packings with continuous size distribution of increasing 
polydispersity. The inset shows that 5(g — >• 0) increases 
continuously with the polydispersity p. This suggests 
that size dispersity is the main factor responsible for the 
numerical results in Rcf. (sj] and the experimental ones 
in Rcf. 0, where p « 0.05 and 5(0) « 0.05. 

For size disperse systems, we consider not only on 5(g) 
but also the compressibility, xt(q)- In matrix notation 
this reads, 



[pk B TxT{q)\ 



x T S" 1 (g)x, 



(3) 



where x T = (x\, . . . , x n ). To find \t one should measure 
all partial structure factors in S(g), and invert this ma- 
trix to get S _1 (g). For a one-component (monodisperse) 
system, one finds p/csTxcKg) = 5(g), and both quan- 
tities are thus fully equivalent. For a binary mixture, 
n = 2, one gets jl4| 



5n(g)5 22 (g)-5 1 2 2 (g) 



x\S-n(q) +x$Sii{q) - 2xix 2 5i 2 (g)' 



(4) 



In Fig. [TJt we show the compressibility measured exper- 
imentally. A clear linear behaviour of the compressibil- 
ity is obtained for low wavevectors. To our knowledge, 
this anomalous behaviour has not been observed experi- 
mentally before. In Fig. we show the compressibility 




FIG. 1: (a) Structure factor and compressibility for the two- 
dimensional packings of disks obtained experimentally. While 
S(q) resembles that of a binary fluid mixture, the compress- 
ibility displays anomalous low-g linear decay, (b) Compress- 
ibility for numerically generated jammed packings of 50:50 
binary mixtures with various size ratios R. All systems share 
the same low-g linear behaviour of the compressibility. 



obtained numerically, using Eq. for binary mixtures 
with different size ratio. For all systems considered, a lin- 
ear behaviour of the compressibility is obtained for low 
wavevectors. This set of results suggests that a relevant 
generalization of Eq. ([!} for binary mixtures is obtained 
by studying xt{q), rather than 5(g). 

While straightforward for discrete mixtures with a 
small number n of components, where S is an n x n 
matrix, the matrix inversion in Eq. (j3J) is conceptually 
and computationally difficult for continuous size distri- 
butions where the size of the S-matrix formally becomes 
infinite. The system studied experimentally in [Io| is of 
this type. To analyze such packings, we derive a system- 
atic approximation for xt{q)- The idea is that if the size 
distribution is sufficiently narrow, xt(q) can be obtained 
perturbatively. To this end, we define ej = (<Ji—cr)/a and 
derive an expansion of xt(q) in powers of the £j up to 
some fixed but otherwise arbitrary order a. To find a suit- 
able starting point for this expansion, recall the relation 
= + pXiXjhij (g) between the partial structure 
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FIG. 2: (a) Zeroth and (b) second order estimate of the com- 
pressibility, Eq. 10, for packings with continuous size distri- 
bution of polydispersity p. While the low-q behaviour of the 
structure factor for a — is featureless (see inset), the com- 
pressibility displays an anomalous linear behaviour at low-g, 
as seen for a — 2. The second order estimate of the compress- 
ibility can be used up to p ~ 30%, while the first order one 
(data not shown) becomes unreliable beyond p ~ 10%. 



factors and hij (g) , the Fourier transforms of the pair cor- 
relation functions gij (r) — 1 @ . The concentration factors 
Xi vary rapidly with ej for narrow distributions, so we 



cannot expand Sij(q). But the pair correlation functions 
gij(r) and hence the hij(q) depend smoothly on particle 
size, so we expand the latter up to e°e°, takinge.g. for 
= 1, hij{q) = h {q) + h 1 (q)(e l + e j ) + h 2 (q)e l e j pf. This 
expansion is inserted into Eq. ([3]) for $y (g) , giving after 
some algebra the compact form 



[pksTxPiq)}- 1 =mJS- 1 (q)m a , 



(5) 



elements ^(g) = 
{0, ...,a}, and cap- 



where the matrix S Q (g) has 
i(e"(q)e"(-q)}, with fi,u 6 
tures fluctuations of the moment density fields e' i (q) = 
Sr=i e i Pifa)- Hence, e°(q) = p(q) is the number den- 
sity field, and S 00 (q) = S(q), the total structure factor. 
The vector = (<5 , ■ ■ • , S a ) in Eq. ([5| has components 
given by the moments of e.; averaged over the particle size 
distribution: 5^ = Y^i=i x i £ i' so that Sq = 1, Si = 0, and 
62 is directly related to the polydispersity, 82 = p 2 ■ The 
result ([5|) relates the compressibility to the S M "(<?) up to 
order a. For a not too large it is simple to compute as it 
only requires the measurement of (a + l)(a+2)/2 reduced 
structure factors. It can be applied to arbitrary particle 
size distributions and is exact for discrete n-componcnt 
mixtures if we choose a = n — 1, as can be shown by 
direct calculation from Eq. ©. 

We have tested our general formula ([5]) using computer 
simulations. When a = 0, one has pk B Tx {0) = S(q), 
which is only exact for n = 1 (monodisperse systems), 
as discussed above. At first order, a = 1, we need to 
invert a (2 x 2) matrix to get ph^Tx? (<?)■ When ap- 
plied to the case of a continuous size distribution this 
formula produces the expected linear behaviour at low q 
for p < 0.10, but deviations appear at larger p. To check 
whether these deviations arc physical, or a result of our 
approximation, we go to second order, a = 2, where the 
required inversion of a (3 x 3) matrix gives: 



J 



pkBTxPiq) 



S 00 S n S 22 + 2S 02 S 01 S 12 - S 00 [S 12 } 2 ~ [S 01 ] 2 S 22 - [S 02 } 2 S n 



S U S 22 _ [5 



I 2 + 2S 2 (S 01 S 12 - 5 02 5 n ) + 5 2 {S m S u - [S 01 ] 2 ) ' 
I 



(G) 



We now find, see Fig. [2}d, that a linear g-dependence is 
obtained for polydispcrsities as large as 30%. This sug- 
gests that the same behavior should in fact be obtained 
for arbitrary size distributions, although measuring the 
compressibility is more difficult when p is very large be- 
cause we need to go to even higher orders a. 

We now discuss the physical significance of our results. 
It is perhaps not surprising, with hindsight, that jammed 
sphere packings have vanishing compressibility since this 
is precisely how the jamming transition was described in 



the opening lines of the paper. However, the quantity 
we call "compressibility" in this work is in fact a par- 
ticular combination of density fluctuations, Eq. (j3|), that 
only reduces to the compressibility at thermal equilib- 
rium when the fluctuation-dissipation theorem holds 0] . 
Remarkably, our results suggest that a similar connec- 
tion between response and fluctuations may exist far from 
equilibrium near jamming, even in real granular packings. 
An obvious connection between response and fluctuations 
holds at q = 0, since xt(0) = (j)~ 1 d4>/dP indeed vanishes 
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near jamming where P <~ (<f> c — (f))^ 1 , so that both sides 
of Eq. <|3j> vanish. It would be interesting to extend these 
considerations to finite q near the jamming transition. 

Our results also illustrate that, for size disperse sys- 
tems, the limit S(q — > 0) is in general not directly re- 
lated to the isothermal compressibility, xt(0) (l5j ]. While 
the latter vanishes in jammed packings, the former is 
free to take any positive value. This explains why previ- 
ous work on size disperse packings failed to observe any 
anomalous behaviour 0, [13] • Both quantities are related 
to the amplitude of fluctuations of the number density, 
but S(q — » 0) ~ (AJV 2 ) captures the total fluctuation 
of N while xt(q — ► 0) quantifies the fluctuations of N 
at fixed composition. This can be seen by defining the 
composition fluctuation fields Cj(q) = Pt(q) — Xip(q) for 
i = 1,. . . ,n — 1. The compressibility from Eq. ([3]) can 
then be rewritten as 

P k B T XT (q) = S(q) - s&^S- 1 (g)soc(g) (7) 

where the (n — l)-dimensional vector so c gathers the 
correlations between number and composition fluctua- 
tions, and the matrix S cc the correlations among the lat- 
ter [l4lri6|. For the compressibility to vanish at jamming 
the two terms must cancel, which implies that local fluc- 
tuations in ./V become fully correlated with composition 
fluctuations. On the other hand, S(q — > 0) remains posi- 
tive because a local fluctuation of A" can be induced by a 
local fluctuation of the mixture composition: fluctuations 
of TV do occur in jammed size disperse packings. 

The behaviour of S(q) in Fig. [2] is clearly inconsistent 
with Eq. (ffj , as noticed previously 0, U3] ■ Does this im- 
ply that jammed size disperse packings are not hyperuni- 
form, as suggested in For point particles, hyperuni- 
formity refers to vanishing bulk fluctuations of the num- 
ber density, as described in the introduction. However, 
for an assembly of spherical particles hypcruniformity re- 
quires vanishing bulk fluctuations of the local volume 
fraction [T3] - While TV and <p are directly proportional 
for monodisperse spheres like those studied in Ref. 
they are not when the packing is size disperse, and thus 
no conclusion can be drawn from S(q) alone. 

A connection between the anomalous compressibility 
studied in this work and suppressed fluctuations of the 
volume fraction in hyperuniform packings can be estab- 
lished. ^From Eq. ©, we realize that whenever the S(q) 
matrix possesses at least one eigenvalue that goes to zero 
at low q, the compressibility vanishes. We have diago- 
nalized S(q) or S a (q) as obtained in our numerical and 
experimental packings, and indeed found that in each 
case anomalous behaviour of the compressibility origi- 
nates from a single vanishing eigenvalue. This implies 
that there exists a particular linear combination of the 
partial density fields which has no bulk fluctuations. A 
detailed analysis of the corresponding eigenvectors shows 
that they are fully compatible with the local definition 
of the volume fraction, 0(q) = (7r/6) Y^7=i ^iPiil)- This 



identification holds exactly for binary mixtures in our 3d 
simulations and in the 2d experimental packings. It also 
holds true for the packings with continuous size distri- 
butions, to the same order in ej that we analyse for the 
compressibility (a = 2, which is accurate up to polydis- 
persity p = 30%). Indeed, we have checked that direct 
measurements of I(q) = (<f)(q)(j>(— q)) coincide with the 
compressibility shown throughout this article in the low- 
q regime where linear behaviour is observed. Therefore, 
the anomalous behaviour of the compressibility reveals 
the absence of bulk fluctuations of the volume fraction. 
This lack of fluctuations was detected also in the recent, 
independent, study of Ref. [Ill ], using a rather differ- 
ent methodology. We conclude that all our size disperse 
jammed packings arc hyperuniform. 

We have demonstrated anomalous behaviour of the 
compressibility in jammed size disperse packings of 
spheres both in simulations and in a granular experi- 
ment, using for the case of continuous size distributions 
an efficient perturbative approach. We have related this 
to suppressed bulk fluctuations of the volume fraction, 
or hyperuniformity (see also [ill ] , thus revealing a struc- 
tural signature of jamming not seen in the conventional 
structure factor. Our work also raises intriguing ques- 
tions about the applicability of fluctuation-dissipation re- 
lations to jammed systems. 
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